Analysis of power-law exponents by maximum-likelihood maps 
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Maximum-likelihood exponent maps have been studied as a technique to increase the under- 
standing and improve the fit of power-law exponents to experimental and numerical simulation 
data, especially when they exhibit both upper and lower cut-offs. The use of the technique is tested 
by analysing seismological data, acoustic emission data and avalanches in numerical simulations of 
the 3D-Random Field Ising model. In the different examples we discuss the nature of the deviations 
observed in the exponent maps and some relevant conclusions are drawn for the physics behind each 
phenomenon. 
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INTRODUCTION 



For the last few decades the study of critical phenom- 
ena has received a great deal of attention in many dif- 
ferent areas of physical'. Criticality is often identified by 
the presence of statistical scale- free distributions of differ- 
ent magnitudes when a system evolves in time or when 
it is driven by an external force. The physical nature 
of the measured response can be very different: energy, 
displacement, magnetization, volume, polarization, resis- 
tivity, etc. In all cases, when the sudden changes (often 
called avalanches) of such magnitudes exhibit a statisti- 
cal distribution compatible with a power law probability 
density function g{x)dx ~ x~°'dx one gains confidence 
about the existence of criticality. In some cases crit- 
icality has also been referred to as "crackling noise" ^. 
At criticality the response of the system is characterized 
solely by the critical exponents a. The theoretical under- 
standing provided by the use of Renormalization Group 
techniques"^ states, in many cases, that critical exponents 
show a certain degree of universality. Thus, it has be- 
come extremely important to determine the values of a 
to a high degree of accuracy and confidence. 

Within this framework, it is important to develop tools 
to test power-law behaviour of data samples and to fit 
critical exponents. The development of these statisti- 
cal tools cannot be done without taking into account 
the limitations inherent to data acquisition. Typically in 
physics, data comes from experiment or computer simula- 
tion and in both cases one has boundaries to the proposed 
scale-free behaviour. These boundaries are not necessar- 
ily sharp or well defined. In experiments one finds un- 
avoidable noise deforming the power-law distribution in 
the small-event region and different kinds of instrument 
saturation in the large-event region. One should take into 
account the fact that it is difficult to find instruments 
(amplifiers, voltmeters, etc.) that allow measurements 
with a range of more than 5 decades. In simulations one 
also finds unavoidable limitations: for instance, one has 
a minimum lattice parameter or particle size that alters 
the small-event distribution and finite-size effects deform- 
ing the large events. Since simulations with more than 



10^ particles are scarce, it is also difficult to find power- 
law distributions extending many decades in numerical 
works. In addition, the existence of deformations in the 
region of small events is understood not only because of 
the reasons discussed above with a physical origin, but 
also because of a mathematical constraint: a pure power- 
law probability density with a > 1 cannot be normalized 
without a theoretical lower limit Xmin- 

It is easy to understand that a naked-eye analysis of 
standard histograms can be easily fooled, not only by the 
lack of statistics (insufficient data in the recorded sam- 
ple), but also by the anomalies in the large- and small- 
event regions. The same may happen with traditional 
fitting methods such as the least-squares method (both 
linear and non-linear) which, for instance, depends on 
the binning process that is performed in order to plot 
the histograms. 

Many years ago most of the community adopted the 
maximum-likelihood (ML) estimation method'*"^ as the 
safest way to treat data, although it is still frequent to 
see papers using alternative, error-prone fitting methods. 
Within this scenario, the work done by M.E.J. Newman 
and co-workers should be pointed out Using ML meth- 
ods they have nicely illustrated how to test the power- 
law character of data and how to obtain good estima- 
tions (and error bars) of critical exponents. One of the 
proposed techniques consists of studying how robust the 
ML exponent is when the analysed data is restricted to 
being higher than an imposed lower cut-off Xiow that is 
varied by several decades. By this method one studies 
the deformation of the fitted exponent due to undesired 
effects in the region of small events. 

In this paper we will study the extension of this tech- 
nique to the analysis of the ML exponent as a function 
of both an imposed lower cut-off Xiow and an imposed 
higher cut-off X^igh- This will render the so-called ML 
exponent maps- . By this method we expect to be able to 
improve exponent estimation in the case in which exper- 
imental or simulation data presents distortions not only 
in the region of small events, but also in the large-event 
region. 

In section |ll] we will revisit the ML method and de- 
fine ML exponent maps. We will include a discussion on 
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numerical methods, evaluation of error bars and analy- 
sis of synthetic data obtained by pseudo-random number 
generation. In section IIIII we will apply the proposed 
analysis technique to the study of three seismological 
catalogues from Japan, the San Andreas fault and the 
very recent activity in the island of El Hierro (Canary 
Islands). In section ITVl we will study three experimental 
cases corresponding to the measurement of the energy of 
acoustic emission (AE) events in different phenomena in 
solids, from previous literature The three sets of experi- 
ments have been carried out with the same experimental 
setup but with different experimental constraints (noise, 
amplification, etc.). The three cases correspond to (i) 
the compression of a porous material (Vycor)^, (ii) a 
cubic-tetragonal structural transition in FePd® and (iii) a 
cubic-monoclinic structural transition in a CuZnAl shape 
memory alloyi^. In Sec. |V]we will illustrate how to ap- 
ply the technique to the study of numerical simulations 
of the 3D Random Field Ising model with metastable dy- 
namics which is one of the prototypical frameworks for 
avalanche criticalit} ^^^'^^ . Finally, in Sec. I VII we will sup- 
ply a summary and draw conclusions. 



II. MAXIMUM-LIKELIHOOD EXPONENT MAP 

Let us consider a sample of measurements {X^} (i ~ 
1, • • • ,N) that we assume to be statistically independent. 
We will denote by capital letters Xmax and X^in the 
smallest and the largest values in the sample set. Our 
aim is to model this set with a power-law probability 
density: 



g{x)dx 



C(7) 



dx 



"\ X ^ X j-i 



(1) 



where 7 is the exponent that we wish to estimate and 
C(7) denotes a normalization function. This normaliza- 
tion function will depend on 7 and the theoretical upper 
(xmax) and lower {xmin) bounds, as will be discussed 
below. (The case Xmax — > 00 is a particular case of 
our analysis). Note that these theoretical limits will not 
necessarily coincide with the maximum Xy^idx and min- 
imum Xmin values in the sample. Nevertheless, given 
the power-law character of the probability density we ex- 
pect that Xmin will be much closer to Xmin than Xmax 
to Xmax- Table U] clarifies the generic definitions of the 
limits and cut-offs that will be used within this section. 
The Likelihood function C is defined as the probability 
that the set of measurements {X^} can be obtained by 
the proposed model: 



N 



(2) 



symbols 


meaning 


^min , ^max 
-^min , Xmax 
Xlow: Xfiigfi 


sharp bounds of the density g{x) 

extreme values in the sample {Xi} 

cut-offs imposed on the sample for the analysis 



TABLE I. Terminology used in this work to define the generic 
bounds and cut-offs used for the theoretical analysis in this 
section 



In order to evaluate the deviations of data with re- 
spect to the proposed model, in this work we will perform 
ML estimations by restricting the original data within 
the imposed lower cut-off Xiow and the imposed higher 
cut-off Xhigh, different from the theoretical limits Xmin 
and Xmax which are generally unknown. We will use the 
symbol 7 to distinguish the exponent estimated within 
a restricted interval from the exponent estimated from 
the whole available sample. We will use n {n < N) to 
denote the number of data of the restricted set. The nor- 
malization factor of the probability density with imposed 
cut-offs is: 



C(7)- 
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high 
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7-1 



(3) 



The best estimation of 7 is consequently found by maxi- 
mizing the likelihood function: 



= 



9 In £(7) 



E 

{X,,^<X,<X^,g^} 



ln(XO-n^ (4) 



where: 



C'(7) _ 1 
C(7) 1 - 7 



high 



X}-7lnX,. 



^high 



^low 



(5) 



We should mention that for the case in which data is dis- 
crete (for instance, in many simulations of lattice models) 
the above treatment should be slightly modified. In this 
case, the data consists of the frequencies of occurrence 
J{k) of a discrete set of values {k} (which we will assume 
to be integers). We would like to fit it with a power-law 
probability function (called Zeta or Zipf)^: 



fc-7 

P{k) — — r kmin l£ k < k„ 

Civ 



(6) 



Following the same procedure as above, when we re- 
strict ourselves to data within imposed cut-offs Kiow and 
Khigh, the normalization function is: 



c(7)= E ^" 



(7) 



k=Ku 



The Maximum Likelihood (ML) estimation method con- 
sists of choosing the value of 7 that maximizes the Likeli- 
hood Function, i.e. the value that makes the sample that 
we have obtained the most likely one to have occurred. 



and the derivative of the Likelihood Function will read: 

d\nC V^'^-'' k-^ ln(k) 

- /(k)ln(k)+jV ^^- ; (8) 
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As opposed to what happens in the case in which one 
considers only a lower cut-off, equations ([5]) and ^ can- 
not be solved analytically. Thus, in this work, we will 
use the false position method in order to find rootsi^. 
This method generates a sequence of recursively smaller 
intervals that always include the root of the equation. 
The monotony of the derivative of In/!^"* ensures that 
the false position method always converges to the root. 
We have chosen arbitrary starting values of 71 = 1 and 
72 = 3.5, and have iterated the algorithm M-times until 
an interval (7M-1 , 7m) is reached with a distance smaller 
than 0.005. 

By changing the Xiow and X^igh cut-offs we can plot 
the values of the ML estimations of 7 using a color 
scale and thus obtain the exponent map. Examples are 
shown in Fig. [T] and throughout the paper. Contour 
lines (in white) will also be shown separating exponent 
values in steps of 0.1. The maps exhibit a triangular 
shape since they are obviously limited by the condition 
Xfiigh > Xiow ■ The main goal of the map is to check the 
existence of a flat plateau (with an homogeneous color 
which is free of contour lines) in which the exponent is in- 
dependent of the cut-offs and thus confirm the scale-free 
behaviour of the data. The map also allows anomalies 
to be identified that can have different origins, as will be 
discussed in the following sections by the use of examples. 

One of the advantages of using the false position 
method for root finding is that we straightforwardly ob- 
tain an estimation of the second derivative of In C at the 



maximum: 



92 ln£ 



9^2 



M 



M-1 



7A/ — J M-1 



(9) 



Assuming Gaussian behaviour of the Likelihood Function 
(which is ensured by the Central Limit Theorem when n 
is large enough) and under very general conditions'^, the 
second derivative provides us with an approximation to 
the standard deviation of the estimated exponent. 
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It is easy to check that this expression, when x 
reduces to the equation proposed in Ref. 0: 

7-1 



(10) 

max ^ CX^i 



(11) 



As a first test of the usefulness of the ML exponent 
maps, we studied synthetic data§. generated according to 
a power-law probability density with theoretical limits 

^min X <^ Xmax ' 



]{x)dx = (7- 1)- 



,1-7 



1-7 

^max 



dx 



(12) 



with 7 = 1.55, Xmin = 1 and Xmax — 10^. The syn- 
thetic data samples were obtained using the RANECU 




10^ lO"* 10^ lo' lo' 10^ 10^ 10^ 
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high 
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FIG. 1. Exponent maps obtained from synthetic data. The 
data sets correspond to a theoretical probability density func- 
tion with exponent 7 = 1.55 and bounds Xmin = 1 and 
Xmax ~ 10*. The label A'^ indicates the size of the sample 
set. The two curved continuous lines indicate the limits above 
which the standard deviation becomes greater than 0.05. The 
red curved line corresponds to the estimation using Eq. (|lip 
and the black curved line to Eq. (|10|) . The vertical black 
straight line marks the highest value in the sample Xmax 



generatorJ^^ for uniform random numbers in the [0,1) in- 
terval and transformed using the method based on the 
inverse cumulative distribution function F~^{z). For 
the proposed probability density the transformation 
function that converts uniform random numbers z into 
the desired ones is given by: 



F-\z) 



z{x. 



1-7 

max 



(13) 



Fig. [T] shows the resulting maps corresponding to four 
samples of increasing sizes N as indicated by the legends. 
The first observation is that the size iV of the sample is 
crucial in order to obtain a clean plateau corresponding 
to the correct exponent. Only when TV > 10^ does the 
plateau extend for several "square" -decades. 

We have indicated the position of the maximum value 
Xmax obtained in the sample by a vertical black line. 
The variations of the fitted exponent observed from this 
line to the right are simply a consequence of changing the 
imposed upper cutoff X^igh in a region without available 
data due to the lack of statistics. This lack of data be- 
tween Xmax and Xfiigh is relevant for the fitting method 
and, a priori, should be taken into account. Furthermore, 
the imposed cutoff X^igh can be moved above the the- 
oretical limit Xmax or set to be oo. This will then be 
equivalent to fitting the synthetic data (generated with 
a theoretical upper limit) with a model without such a 
limit. In the maps in Fig. [T] we have kept the imposed 
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low low 



FIG. 2. Comparison of the behaviour of the exponent fitted 
to synthetic data as a function of the lower cut-off Xiow when 
Xhigh is fixed. Red squares correspond to X^igh = Xmax ~ 
10* (the theoretical limit of the synthetic data), green circles 
to Xhigh = Xmax (the maximum value found in the sample) 
and blue triangles to Xhigh = oo. 

cut-off Xhigh below the theoretical cutoff Xmax , but note 
that for real data (e.g. in the following examples) this 
theoretical cut-off will be unknown. 

Fig.[2]shows the behaviour of the fitted exponent 7 as a 
function of the lower cut-off Xiow for the same synthetic 
data samples as in Fig. [TJ We have compared three dif- 
ferent ML estimations corresponding to three choices of 

Xhigh- 

1. The green circles correspond to the ML estima- 
tion of the exponent obtained by fixing Xhigh — 
Xmax = 10^ , the theoretical upper limit. This is 
nothing more than the profile of the ML exponent 
map along the vertical right border at 10® in Fig.[T] 
This would be the correct way to perform the ML 
estimation if one could have a priori knowledge of 
the true upper cut-off. 

2. The blue triangles (higher symbols) correspond to 
the ML estimation obtained by neglecting the ex- 
istence of an upper limit in the power-law distribu- 
tion, i.e. fixing Xhigh = 00. This would precisely 
correspond to the method proposed by Ref. H. 

3. The red squares (lower symbols) correspond to the 
ML estimation obtained by fixing the higher cut- 
off of the data to the maximum value found in the 
sample Xhigh = X,nax- It corresponds to the profile 
of the map along the vertical black lines in Fig. [TJ 

For small N Method 3 underestimates the exponent 
because it neglects the fact that no data has been ob- 
served between Xmax and Xmax- However, for > 10^, 
one can see that Method 3 renders exponents that are 
very similar to the correct ones (Method 1). For N — 10^ 



Method 3 is clearly better than Method 2, which neglects 
the existence of an upper boundary to the distribution: it 
can be seen that the two coinciding estimation methods 
(red squares and green circles) exhibit a larger plateau 
(by more than 1 decade) than the blue triangles. 

It is also interesting to discuss the strong fluctuations 
of the ML exponent close to the diagonal of the maps, 
which can be observed in Fig. [T] These are due to 
the small size n of the restricted sample set that in- 
creases the statistical fluctuations. We can locate the 
region where the standard deviation of the ML estimate 
is lower than ±0.05^^ by using the formulas pUj) or (ITT|) . 
These low-error regions correspond to the areas below the 
black- and red-curved lines, respectively. As can be seen, 
both estimations of the error differ, even in the region 
of small events. The estimation that takes into account 
the existence of both a low and a high cut-off obtained 
from Eg. ([To)) gives a better separation between the re- 
gions with meandering contour levels and the smooth, 
flat plateau. 

To conclude this section, let us summarize what we 
have learned from the analysis of synthetic data: it does 
not make much sense to increase the higher cut-off Xhigh 
above the maximum value in the data sample Xmax un- 
less we have independent information of the theoretical 
limit Xmax- Therefore, in the maps presented in the fol- 
lowing sections we will scan the higher and lower cut-offs 
only within Xmin and Xmax- Thus, the vertical right 
border of the maps will coincide with the vertical black 
line plotted in Fig. [T] In addition, we will use the error 
estimation proposed by equation ([TO]) and plot, on the 
maps, the curved line separating the region with error 
bars greater than ±0.05 (above the line) from the region 
with lower error bars (below). 



III. SEISMOLOGICAL CATALOGUE ANALYSIS 

The Gutcnberg-Richter lawi^iiS describing the sta- 
tistical distribution of earthquake magnitudes is one 
of the most famous examples of a scale-free phe- 
nomenon already discussed using ML methods in pre- 
vious works^'^. Theoretical studies have proposed 
different physical models (Burridge-KnopofSi, Olami- 
Feder-Christensen^^, damage rheology^^,...) and frame- 
work theories such as the so-called Self Organized 
Criticalitj«2ir— which have explained, to a certain ex- 
tent, the reasons behind this critical behaviour. It is not 
our purpose to gain any understanding of seismology, but 
only to use some of the available earthquake catalogues 
in order to test the behaviour of ML exponent maps. 

Earthquakes are historically characterized by a quan- 
tity called magnitude M, which aims to be a logarithmic 
measure of the "size" or "energy released" during the 
earthquake. The Gutenberg-Richter lawi^ refers to the 
number of earthquakes A'^> (M) with a magnitude larger 
than a certain value M As a function of magnitude, the 
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Gutenberg-Richter law can be written as: 

Ny{M)(x-bM (14) 

with b ^ 1.0. The measurement of earthquake magni- 
tudes and energies is still a challenging issue for seismol- 
ogy. As the medium is too large to collect a significant 
amount of radiated energy, seismologists must rely on 
measurements of the sparse network of seismic stations 
in order to locate and estimate the "size" of an earth- 
quake. Because of this, many different criteria are used, 
depending on the region, earthquake energy range, avail- 
able instruments, etc. The catalogues usually contain 
mixed data corresponding to different definitions of the 
magnitude M. These definitions are not fully equiva- 
lent, especially for small earthquakeai2i22. There are also 
different definitions of the "energy" associated with an 
earthquake (all of them are approximately linearly re- 
lated): seismic moment, strain energy drop, radiated en- 
ergy, etc.. In this work we will use a broadly accepted 
formula^^ that allows an approximate conversion of the 
different "magnitudes" M to the minimum strain energy 
drop E as: 

Iogio£: = 1.5M-I-4.8 (15) 

where E is the energy in Joules. Using the Gutenberg- 
Richter law ((T4)) and Eq. (|15p . one can write the prob- 
ability density for earthquakes with energies between E 
and E + dE as: 

p{E)dE - E-'dE (16) 

where e = 1 + (6/1.5) ~ 1.67 is the expected exponent 
characterizing the power-law distribution of earthquake 
energies. 

We have computed the ML exponent maps correspond- 
ing to three earthquake catalogues: 

1. The subduction process taking place in the Japan 
Trench makes it one of the most active seismo- 
logical regions in the world. The area is quite 
well documented because of the vicinity of the 
Japanese islands. We studied the exponent map 
corresponding to the energy distribution for all 
the seismological events registered as earthquakes 
in the ANSS29 catalogue from 2000/01/01,00:00:00 
to 2011/11/09,17:32:36 within the region enclosed 
between latitudes 28°N and 48°N and longitudes 
128°E and 148°E. The registered data correspond 
to the = 14509 events above M — 2.7, where 
the the Tohoku earthquake of 2011/03/11 was the 
most serious event with an estimated magnitude of 
M = 9.0. 

2. San Andreas fault system^S, beneath the region oc- 
cupied by the states of California and Nevada, there 
is probably the most frequently monitored seismic 
region in the world and the best documented in 
catalogues. We will therefore take it as a precise 



example of a seismological strike-slip process. The 
data analysed here corresponds to the seismic sig- 
nals registered as earthquakes in the ANSS^? cata- 
logue with its epicentre within the area of latitudes 
between 30°N and 42°N and longitudes between 
114°W and 126°W during the period between 
2000/01/01,00:00:00 and 2011/11/09,17:43:00. In 
order to avoid the presence of a possible noise back- 
ground, we selected only those earthquakes with a 
magnitude greater than M — 0.4. The strongest 
earthquake was recorded on 2005/06/15 off the 
Coast of Northern California with a magnitude of 
M = 7.2. The data set has N = 453372 events. 

3. As a completely different seismological phe- 
nomenon (very localized in space and time), we 
considered the recent submarine volcanic eruption 
of La Restinga off the island of El Hierro (Ca- 
nary Islands), which started in summer 2011. The 
volcanic activity triggered an earthquake swarni^ 
which is expected to have quite different behaviour 
from typical tectonic processes. We considered 
the data obtained from the IGN^^ catalogue from 
2011/06/08, 00:52:00 until 2012/2/07, 12:00:00 in 
the region enclosed by latitudes from 26.8°N to 
27.6°N and longitudes from 17.85°W to 18.2°W. 
The data set has N = 12158 events 

The three ML exponent maps are shown in Fig. |31 We 
have kept the same scales in order to clearly reveal the 
different size N of the statistical samples. Fig. |4] shows 
the behaviour of the fitted exponent e as a function of the 
lower cutoff Eiow when the higher cut-off Ehigt is fixed 
to the maximum value in the sample set, i.e. the profile 
of the map along the vertical right borders in Fig. [3l 

The first observation is that there is an almost perfect 
plateau exhibited by the San Andreas data (middle dia- 
gram) for a value close to the expected theoretical value 
e = 1.67. Despite some deformation, indications of a 
plateau are also observed for the other two sets, Japan 
(top diagram) and El Hierro (bottom diagram). For the 
El Hierro data the coincidence is remarkable, given the 
different physical origins of the earthquake sequence. 

A second important observation is the deformation of 
the plateau (towards low exponent values) for the Japan 
data in the region of Eio^ < 10^° J. A plausible ex- 
planation for this deformation is that the statistics for 
small earthquakes in the Japan catalogue is incomplete. 
The same tendency can be observed for the San Andreas 
data, but for much lower minimum cut-offs Eiow < 10* J, 
almost coinciding with the lower border of the map in 
Fig. [3] (middle diagram) . The oscillation of the fitted ex- 
ponent for the Japan data that can be seen on the maps 
(as contour lines with a parabolic horizontal shape start- 
ing from the right border) is also surprising as well as 
the maximum (about Eiow ~ 10^^ — 10^^ J) on the profile 
shown in Fig. |4l We are unable to provide an explana- 
tion for this behaviour, but it could be caused by the 
different methods used to estimate magnitudes and/or 
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FIG. 3. ML exponent map corresponding to the earthquake 
data from Japan, the San Andreas fault and El Hierro. White 
contour lines are separated by 0.1 units. The region above 
the black line corresponds to estimated statistical error bars 
greater than ±0.05. 



energies depending on the earthquake magnitude range. 

IV. ACOUSTIC EMISSION DATA ANALYSIS 

As a second set of experimental examples, we have 
focussed on much smaller energy scales. Different pro- 
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FIG. 4. Behaviour of the ML exponent as a function of the 
lower cut-oiT Eio^u for a fixed higher cut-off equal to the max- 
imum value in the sample set. Only a few error bars are 
indicated. The horizontal dashed lines show the theoretically 
expected value 1.66 and an error bar ±0.10. 



cesses, which have been classified as critical or crackling 
noise, take place in solids, exhibiting a certain degree of 
disorder when driven by an external force or by a tem- 
perature ramp. Examples include superconductivity'^5,, 
capillary condensation'^*, acoustic emission in structural 
.j Q transitions^"—, Barkhausen noise in magnetism^i^, 
fracture^, etc. In many cases, theoretical studies have 
provided general frameworks to understand the origin of 
this criticality^i^^i^. 

The following three experimental examples correspond 
to data recorded using the Acoustic Emission (AE) 
technique*® . Propagation of cracks or the sudden move- 
ment of internal interfaces generate acoustic waves in the 
ultrasonic range that propagate through the solid and 
that can be recorded by appropriate transducers on the 
surface. The AE method is equivalent to the method used 
to monitor earthquakes, but on a much smaller scale. It 
is interesting to describe here some details in order to 
understand the deviations that will be observed in the 
maps. 

The most common piezoelectric transducers, coupled 
to the sample surface, generate a voltage signal V(t) pro- 
portional to the speed of the incident elastic wave that is 
amplified. Using a predefined threshold (above the un- 
avoidable experimental electrical and mechanical noise), 
it is possible to define individual AE events^. The be- 
ginning of an event occurs at a time ti when the voltage 
exceeds the threshold. The end of the event occurs when 
the signal falls below the threshold at t2 and remains 
below the threshold for more than a certain pre-defined 
time called the Hit Definition Time (HDT), typically in 
the range of 10-100 fis. The fast integration of the V^{t) 
signal from ti to t2, normalized by a reference resistance, 
renders an estimation of the energy recorded by the trans- 
ducer which is assumed to be proportional to the energy 
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released by the physical process generating the elastic 
wave. 

It is worth mentioning some experimental limitations 
of most standard set-ups: (i) acquisition systems, due 
to limited memory, have an internal maximum limit on 
the duration of a signal as well as a maximum limit on 
the voltage that saturates the amplifier. In the case that 
such maxima are exceeded, the signal is truncated both in 
voltage and/or duration. This represents a deformation 
in the large-event region and a not totally sharp cut-off in 
the measured energy since both voltage and duration can 
be independently exceeded, (ii) A second experimental 
problem that needs to be considered, is due to the at- 
tenuation of ultrasounds inside the material and the dis- 
tance from the source of the AE event to the transducer. 
If the studied samples are small (compared to typical 
length scales for exponential attenuation or compared to 
typical transducer sizes) we expect that data recorded 
by a single transducer would not be very distorted by 
the distance to the source. However, if samples become 
large, the quality of the overall power-law distribution 
becomes poorer. One could then use several transduc- 
ers to locate the position of the source of the event and 
correct for attenuation. This has been achieved in some 
cases^i^, thought, in general, it is a complicated pro- 
cedure. The examples below correspond to the use of a 
single transducer. (iii)A third problem is that counting of 
the small signals is lower than expected for different rea- 
sons: an important fraction of signals that happen to be 
very short in time and/or amplitude cannot be detected 
by the acquisition set-up (they are too short for the sam- 
pling frequency or the amplitude is below the threshold 
value), or some of the small signals may be overlapped 
by the tails of previous signals due to dead-time HDT. 

In the cases analysed in the following three subsec- 
tions, the same experimental set-up for the acquisition of 
AE has been used: a PCI2-system from the MISTRAS 
Group, which consists of an 18-bit A/D converter work- 
ing at a base sampling rate of 40 MHz. The transducers 
are also the same in the three examples (microSO). This 
makes the ML exponent maps easy to compare. The 
amplification factor, the threshold and the number of 
recorded signals used in each case are different, given 
the different noise conditions and differing nature of the 
studied phenomenon and driving force. Although the ac- 
quisition card is only 18 bits (~ 5 decades), the measured 
energies by the fast integration algorithm may theoret- 
ically extend many more decades. In order to compare 
the different studies, we will keep the same scales: from 
10^^^ J to 10^^ J for energy and from 1.0 to 2.5 for the 
fitted exponent e. 



A. Mesoporous Si02 under compression 

The study of noise in porous materials under compres- 
sion is important for the prediction of accidents in min- 
ing. In this first example* two samples of a porous mate- 



rial Si02 (Vycor) with two parallel faces are compressed 
between two plates with a lineally increasing load in time. 
The AE sensor is attached to one of the plates. There are 
two sets of data studied corresponding to the AE events 
recorded at two loading rates: 0.2 kPa/s and 1.6 kPa/s. 
The study of the infiuence of the driving rate is impor- 
tant because, in some cases, it affects the fitted exponents 
due to the overlap of large and small avalanches^i^. The 
material cracks under compression and the recorded AE 
signals show a power-law distribution of energies. In this 
case a preamplifier of 60 dB was used and the threshold 
was selected to be 26dB. The number of recorded signals 
was N = 11022 and N = 28652 for the first and second 
set, respectively. 

Fig. [S] shows the ML exponent maps obtained with the 
two driving rates. It is clear that, in both cases, there is 
a vast region of the map with a constant plateau corre- 
sponding to a common value of the exponent close to 1.4. 
This means that the exponent is robust against changes 
in the cut-offs by several decades. The small coloured 
spots observed close to the diagonal of the map as well 
as the large peak in the upper right-hand corner cor- 
respond to expected statistical fiuctuations, above the 
black line that indicates when the statistical error bar 
becomes larger than ±0.05. The black region in the bot- 
tom left-hand corner corresponds to the fact that the 
exponent decreases when only an important fraction of 
very low signals is included in the analysis which gives 
an erroneous estimation for the reasons explained above. 
Basically, no other deviations are observed in this exam- 
ple. 

Fig. [S] shows the analysis of the profiles of the map 
along the right vertical border (Ehigh = Emax) compared 
to the analysis proposed in Ref. la assuming no upper 
limit for the power-law distribution {Ehigh = oo). As 
can be observed, the two profiles are similar, but are 
not identical. They show a clean plateau at 1.39 ± 0.02 
for three decades (10~^^ - 10~^^). The profile obtained 
by fixing the highest cut-off to the maximum value in 
the sample (i.e. along the right-hand border of the ML 
exponent map in Fig. [5]) shows fiatter behaviour along, 
at least, one more decade. 

Therefore, one can conclude that this example shows 
a very robust power-law behaviour, comparable to that 
observed for seismological data. 



B. Structural transition in FePd 

In this second case we present an analysis of the AE 
recorded during a martensitic transition in Fe-Pd2. This 
alloy has received a lot of theoretical interest due to the 
simple symmetry relation between the high-temperature 
phase (cubic) and the low-temperature phase (tetrag- 
onal). The transition can be induced by applying an 
external stress or by changing temperature. Two sam- 
ples were considered in the study: a single crystal and a 
polycrystal with the same composition Fe68.8Pd31.2- The 
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FIG. 5. ML exponent maps corresponding to the AE recorded 
during the compression of Vycor samples at two different 
rates, 1.6kPa/s (above) and 0.2 kPa/s (below). White lines 
correspond to the exponent contour levels with a separation 
of 0.1. The region above the black line indicates an estimated 
error bar grater than ±0.05. 



transition was induced by changing the temperature at 
different driving rates between O.lK/min and lOK/min. 
Apart from the critical distribution of energies during 
the events, the AE study revealed other interesting fea- 
tures (which were not found by other techniques). The 
transition, on cooling, started at the same temperature 
Ms = 2A6±2K for the two samples. As soon as the tran- 
sition started the formation of new tetragonal domains or 
the advance of interfaces separating previously nucleated 
domains, generated AE events. Due to the thermoelastic 
behaviouriS of the transition the AE extended for many 
degrees (~ 35K) until the sample was fully transformed. 
On heating, the reverse process occurred with low ther- 
mal hysteresis (~ IK). 

In this case AE was amplified 60dB and the threshold 
was chosen to be 22dB. The four sets of data that we 
will analyze here correspond to a driving rate of IK/min 
for both the single crystal and the polycristal and for 
cooling and heating runs. In order to increase statistics. 
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FIG. 6. Profiles along the right vertical border of the exponent 
maps corresponding to the Vycor experiments, compared to 
the profiles proposed in Ref. [y. Some typical estimated error 
bars are shown. The horizontal line indicates the value e = 
1.39 proposed in Ref. §. 



data were accumulated over 20 ramps. This accumu- 
lation technique would correspond to a real increase of 
the statistics only if the data recorded on each ramp is 
independent of the previous data. This assumption is 
doubtful in the case of structural transitions since it has 
been demonstrated that samples exhibit a learning pro- 
cess that increases the correlation of the signals between 
consecutive ramps^. In this case the ramps correspond- 
ing to the same driving rate were not strictly consecutive 
since other driving rates were added in between, thus the 
independence of the recorded data is not clear. 

The total number of recorded signals was A'' = 171056 
(single crystal, cooling), N = 111840 (single crystal 
heating), N — 192596 (polycrystal cooling) and N = 
58501 (polycrystal heating). Typically, the number of 
AE events is not symmetric during heating and cool- 
ing ramps and there is a lack of understanding of this 
phenomenor*^. 

The studyS. concluded that in the four cases the en- 
ergies of the individual AE events were power-law dis- 
tributed ( ^ 4 decades for the single crystal and ^ 3 
decades for the polycrystal). The exponents were fit- 
ted using the ML method and were almost independent 
of the heating/cooling rate. The values of the cut-offs 
were selected in the region where the statistics were suf- 



9 



ficiently high. For the single crystal both cooling and 
heating ramps exhibited an exponent compatible with a 
value of 1.64 ±0.10. For the polycrystal a clear deviation 
of the exponent was found when the heating (2.0 ± 0.1) 
and cooling ramp (1.59 ± 0.10) were found. So far, we 
have no explanation for this deviation. 

Fig. [7] and [8] show the ML exponent maps correspond- 
ing to the four cases. The graphs above correspond to 
cooling data and the graphs below to heating data. As 
can be observed the coloured spots attached to the diago- 
nal, which correspond to statistical fluctuations are larger 
(in absolute terms) than in the previous example. This 
suggests that, although the recorded number of signals is 

10 times larger, most probably data corresponding to 
the different 20 runs were correlated and did not effec- 
tively increase the statistics. Thus, the estimated error 
bars were, most probably, underestimated by a factor of 
-\/20. Therefore, in order to plot the curved black line 
on the map separating the zone with large error bars, we 



have required that y/20aj = 0.05. It can be seen that, 
indeed, the curved line separates the zone with fluctuat- 
ing contour lines from the flat region. By observing the 
map corresponding to the heating runs (Fig. |S] bottom) 
for the polycrystalline sample one clearly sees that the 
large fitted exponent ^ 2.0 is doubtful and most proba- 
bly is a consequence of the lack of statistics. The region 
with small enough error bars (below the black line) is 
very small and it is difficult to identify any plateau. For 
the cooling ramps (Fig. [S] top) the plateau is not very 
clear but at least the separation of the contour lines is 
much wider. 

Fig. |9] shows the corresponding profiles along the right 
vertical border of the map. Error bars have been in- 
creased by a factor of V20 compared to the ones reported 
in the original paper—. The dashed lines indicate the pro- 
posed exponent 1.64±0.1. The hypothesis that this might 
be a common value for the exponent corresponding to the 
four cases cannot be ruled out due to the lack of statis- 
tics. More measurements for the polycrystalline sample 
would be required to fully clarify this point. 
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FIG. 7. Exponent maps for AE in the FePd single crystal 
corresponding to heating ramps (below) and cooling ramps 
(above). The region above the black line corresponds to es- 
timated error bars greater than ±0.05, after correcting by a 
factor of y/20 due to the possible correlations between mea- 
surements (see text). The contour lines (in white) are sepa- 
rated by 0.1 units. 



C. Structural transition in CuZnAl 

This third example corresponds to a recent study of 
AE during the martensitic transition in a CuZnAl shape- 
memory alloy In this case the transition is also ther- 
mally induced from a monoclinic multivariant structure 
at low temperatures to a cubic structure at high temper- 
atures. The purpose of the analysis performed in Ref. 
[Tol was not to demonstrate the power-law distribution 
of the AE events (which had been shown previously in 
different studiesS^'^^'^), but to compare the observed expo- 
nent with the exponent corresponding to the very large 
avalanches that can also be recorded by a sophisticated 
calorimetric study. The sample studied was, therefore, 
larger than those studied previously. This, as explained 
above may lead to distortion of the power law. Further- 



more, a larger sample typically requires a more powerful 
temperature control set-up, involving higher electric cur- 
rents and thus involving greater noise. Since the study 
was focussed on large avalanches and noise was high, the 
amplification was set to a much lower factor, 40 dB, and 
the threshold to 45 dB. The data set analyzed consists 
of iV = 17936 signals corresponding to a unique heating 
ramp. 

As can be seen in Fig. llOi the exponent map shows 
a poorer quality without any clear plateau. The pro- 
file shown in Fig. Ill) corresponding to the right vertical 
border of the map, suggests an incipient plateau around 
e — 2 for slightly more than one decade, but this is not 
fully conclusive. Thus, the exponent map analysis in this 
case indicates that either the number of signals is too 
small to draw conclusions about the power-law behaviour 
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FIG. 8. Exponent maps for AE in the FePd polycrystal 
corresponding to heating ramps (below) and cooling ramps 
(above). The region above the black line corresponds to es- 
timated error bars above 0.05, after correcting by a factor of 
y/20 due to the possible correlations between measurements 
(see text). The contour lines (in white) are separated by 0.1 
units. 
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FIG. 9. Exponent values versus lower cut-off Eiow for AE in 
FePd. The horizontal lines show the value 1.64±0.10 proposed 
in Ref.- 
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FIG. 10. ML Exponent map corresponding to the martensitic 
transition in a CuZnAl sample during a heating ramp. 



of the data or the distribution shows a deformation prob- 
ably due to the fact that the sample was too large and 
attenuation introduces a length scale. Note, however, 
that we are not fully compromising the results pointed 
out Ref. [13 • The conclusions were based not only on 
AE signals, but also on the coincidence of the observed 
incipient plateau at e = 2 with the plateau observed by 
a different experimental technique. 



V. SIMULATION DATA 

As a last case we will analyze data corresponding to 
numerical simulation of a lattice model. This is an in- 
teresting illustrative case of the advantages of ML maps 
compared to the previous proposed analysis assuming 



no upper limit. This is because the distortions affect- 
ing the large avalanche region are not due to measure- 
ment problems, but are intrinsic to the finite size of 
the model. The study corresponds to the 3D Gaus- 
sian Random Field Ising Model (3D-GRFIM) driven by 
an external field H(t) with metastable dynamics at zero 
temperature^^. The model is based on the original Ising 
model with the addition of random internal fields hi act- 
ing on each spin 5*^ = ±1. The values hi are quenched 
and distributed according to a Gaussian probability den- 

sity p{h)dh = --^=j^e^^dh with zero mean and variance 

R^. The parameter R is usually referred to as the amount 
of disorder in the system. The Hamiltonian of the system 
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FIG. 11. ML fitted exponent as a function of the lower cutoff 
Eiow when considering a higher cut-off equal to the maxi- 
mum recorded signal in the data sample. The horizontal lines 
indicate the value 2.15 ± 0.10. 



reads: 

N N 

H = -^5^S',-^/i,5^-i?(i)^5, (17) 

where the Si spin variables are defined on a regular cu- 
bic lattice and the first sum extends over all nearest- 
neighbour pairs. The simulations of the model start from 
a saturated configuration {Si} = —1 and H = — oo . The 
field is then adiabatically increased and the spins flip ac- 
cording to the local relaxation rule: 

S^ = sign Sj + h, + H{t)j (18) 

where the sum extends over all the z = 6 neighbours 
of the spin Si . With this metastable dynamics the sys- 
tem evolves following a sequence of magnetization jumps 
(avalanches) occurring at certain fixed values of the ex- 
ternal field separated by periods of inactivity in which 
the field is increased without producing any spin flip. 

The model has been widely used as a prototype model 
for the study of avalanche dynamics. It has been success- 
ful in explaining different features of the magnetization 
process in ferromagnets: the presence of rate indepen- 
dent hysteresis, the return point memory property and 
the existence of Barkhausen noise^^. Extensions of the 
model have been also used for the understanding of other 
athermal first-order phase transitionSi ^^'^^ 

Here we will focus our attention on the distribu- 
tion D{s,R) of the sizes s (number of flipped spins) of 
the avalanches obtained along the magnetization process 
from H = —oo to H = +oo, i.e. the so-called integrated 
distribution. In the thermodynamic limit and when the 
amount of disorder R is tuned to the critical value Rc, 
this distribution is expected to be a power la w^^'^^~ — 



characterized by a critical exponent called t' = t + a/3S 

D{s,R,) = ' (19) 

For values of disorder above i?c it is expected that the dis- 
tribution of avalanche sizes is exponentially damped and 
thus, in the thermodynamic limit, the discontinuities in 
the magnetization Am — s/L^ vanish. Below Rc, it is 
expected that the distribution of avalanches is also expo- 
nentially damped, but there will exist a unique massive 
avalanche with a size proportional to that is respon- 
sible for a magnetization discontinuity as should occur 
given the first-order character of the phase transition. 

In the numerical simulations on a finite lattice (Lx Lx 
L) with periodic boundary conditions, nevertheless, the 
distribution of avalanche sizes D{s,R,L) behaves quite 
differently. Several effects deform the power-law charac- 
ter at the "pseudo" -critical point: 

On the one hand, avalanche sizes are limited from 
above by the finite size s < L^. The fact that in all 
critical phenomena with an associated diverging correla- 
tion length, distortions occur well below the limit is 
well known. Among other reasons, close to the critical 
point Rc, the avalanches are expected to be fracta l "'^^1^° 
and thus exceed the lattice side L when its size is much 
smaller than L^. Such avalanches that expand the lat- 
tice side L in, at least, one dimension, are the so-called 
spanning avalanches. For a cubic lattice with periodic 
boundary conditions, they can be easily classified as 
ID-spanning, 2D-spanning or SD-spanning depending on 
whether they span the lattice side in 1,2 or 3 spatial di- 
mensions. 

On the other hand, some small avalanches (sometimes 
called lattice animals within percolation theory^) may 
occur for probabilistic reasons even very far away from 
the critical point. For any value of R and H it is possible 
to compute the probability that two neighbouring spins 
(surrounded by negative spins) flip simultaneously. The 
same computation can be carried out for the probability 
of a group of three, four, etc. neighbouring spins with 
a certain topological configuration to flip simultaneously. 
It is clear that in such a computation the number of con- 
figurations in which spin clusters (lattice animals) with 
a certain size s can be constructed plays a fundamental 
role. Such a phenomenon has (a priori) no connection 
with the critical point. The total number of such non- 
critical avalanches is, nevertheless very large (it increases 
with L^) and renders a distribution of small avalanches 
that has general, fast decreasing behaviour with s, but 
is non-monotonous. At the critical point this effect over- 
laps with the proposed power-law distribution and thus 
shows a deformation in the small-size region of the distri- 
bution D{s; R, L), which typically can be observed below 
s ^ 10. Our goal here is to show how these distortions at 
large and small sizes show up on the ML exponent map. 

Before the discussion of our results, it is important 
to recall that the exact values of Rc and t' have not 
been definitely established. The problems are precisely 
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due to the fact that it is difficult to deal with spanning 
avalanches in the simulations on a finite lattice. There 
have been two main approaches: 

1. Dahmen and coworkers found i?^ = 2.16 ±0.03 and 
r' = 2.03 ± 0.03 by performing a scaling collapse of 
the avalanche distribution, neglecting the fact that 
there should be a dependence on L in the distribu- 
tion. This was first done^ on very large systems 
(up to 1000'^) but averages over very few realiza- 
tions of disorder. The collapses clearly revealed 
that corrections to scaling were needed. Later^- 
similar collapses were done which included a unique 
system size 64^. In both cases, the data were re- 
stricted to amounts of disorder 2.25 < i? < 6.0 
well above R^. By this method they avoided most 
spanning avalanches in the simulations but paid the 
price of working too far from the critical point and 
so they had to extrapolate the value for the expo- 
nent. 

2. Furthermore, there have been studies precisely fo- 
cussed on the behaviour of spanning avalanches 
which analyse how they concentrate close to the 
critical poin t^^'^° . Such studies have performed 
finite-size scaling analysis of the number of span- 
ning avalanches and have obtained a higher value 
of the critical amount of disorder _Rc = 2.21. To 
do so, they proposed a method (called method-2 in 
Ref. 160 ) to separate the 3D-spanning avalanches 
that will correspond to the massive avalanches 
that clearly disturb the power-law distribution of 
avalanche sizes in the region of large events, close 
to the critical point. Such 3D-spanning avalanches 
are identified when they are the unique spanning 
avalanches in the full-field excursion. This allows 
them to be filtered from the statistical analysis. 
The reason behind the separation method is that 
when massive avalanches occur, they fill such a 
large fraction of the system that they do not allow 
for any other spanning avalanche to take place. It 
was shown that this filtering method, although not 
perfect, gave more consistent results in the finite- 
size scaling analysis than the method of discarding 
all the spanning avalanches. 

The differences between the two above approaches are 
even more subtle and difficult to summarize here. Let us 
simply remark that the two analysis use a slightly differ- 
ent scaling variable to measure the distance to the critical 
point and that the exponent r' — 2.03 ±0.03 proposed in 
Ref. i58i is interpreted not as a true critical exponent, but 
as an effective exponent r^^j in Ref. [H. Instead, they 
propose that close to the critical point the distribution 
of avalanche sizes (neglecting the massive avalanches) 
will be dominated by non-critical lattice animals (in the 
small s region) and by the so-called non-spanning critical 
avalanches (in the intermediate and large s region). Only 
the last kind of avalanches may exhibit true power-law 
behaviour with an exponent r„sc = 1-65 at Rc = 2.21. 



We have performed numerical simulations with system 
sizes ranging from L = 32 to L = 256 and values of R 
within the range 2.15 < i? < 2.24. For every size L and 
every i?, averages were taken over 2000 configurations of 
the random fields for L < 128 and over 1400 for L = 256. 
We used the RANECU random number generator and 
a Box MuUer-Polar Marsaglia algorithm to generate the 
Gaussian random fields. Simulations of the metastable 
dynamics were done with the sorted list algorithm^. The 
typical sizes of the sample sets of recorded avalanches 
(corresponding to each single value of R) range from 10^ 
for L = 32 to 10^ for L = 256. 
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FIG. 12. ML exponent maps corresponding to the 3D-GRFIM 
for R = 2.20 and increasing system sizes from L = 64(below) 
to L = 256 above. The left-hand column shows the maps that 
where obtained by considering all the recorded avalanches, 
whereas the right-hand column, the method proposed in 
Ref. i60| was used to suppress the massive, large non-critical 
avalanches. The region above the black line (very close to the 
diagonal of the map) corresponds to an estimated error bar 
of ±0.05. Contour lines in white are separated by 0.1 units 

Fig. [T^ shows the exponent maps corresponding to 
R = 2.20 and increasing values of L. The first column 
corresponds to the exponent obtained by considering all 
the avalanches and the second column to the analysis 
after the suppression of the massive avalanches, as ex- 
plained in Ref. 60. As can be observed in the first col- 
umn the maps display an approximately triangular black 
region in the upper right-hand corner that also extends 
along the right edge which is due to including of such 
"massive" avalanches that they strongly distort the pos- 
sible power-law distribution. Massive avalanches over- 
populate the large s region and thus decrease the value 
of the fitted exponent. Such a deformation disappears 
when the data is filtered, as can be seen in the right col- 
umn in Fig. 1121 
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FIG. 13. Profiles of the ML exponent maps along the vertical 
right border for different values of R and L. Data in the left 
column corresponds to the full set of avalanches. On the right 
column data corresponds to the same sets after filtering the 
massive avalanches following the method proposed in Ref. H^. 



Close to the bottom horizontal boundary of the maps, 
exponent oscillations (seen as a sequence of parallel 
horizontal contour lines) can also be observed in all 
the graphs. These oscillations are due to non-critical 
avalanches (lattice animals) which are always present and 
difficult to subtract. 

Besides these two observed deformation regions, it is 
not clear that a clean plateau exists close to the bottom 
right-hand corner, with an area increasing with L. The 
maps show a slow decrease of the exponent from 2.1 to- 



wards 1.7 when increasing the lower cutoff Siow, and this 
effect does not seem to disappear for large system sizes. 
Instead, a region with a value of the exponent close to 
1.6— 1.7 seems to develop close to the upper corner when 
the system size is increased up to i = 256. This effect 
is clearer when the data has been filtered (right column) 
and suggests that a true power law may only be observed 
when small avalanches are not included in the analysis. 

Fig. [13] shows the profiles of the maps obtained by fix- 
ing Sfiigh — Smax, for different values of R between 2.16 
and 2.25 and L from 32 to 256. The column on the left 
corresponds to the analysis of all the avalanches and the 
column on the right to the analysis after filtering the 
massive avalanches. Again, it is difficult to find any clear 
evidence of a plateau growing with L. The only exception 
seems to be for R = 2.20 — 2.21 in the region of inter- 
mediate and large cut-offs (10** < siow < 10^), where a 
plateau seems to form and become broader when L in- 
creases. This plateau has a height approaching r' ~ 1.7. 
It can be observed not only on the filtered data (right 
column), but also as an inflection when the avalanches 
are not filtered (left column). 

The profiles in Fig. [T51 also allow the oscillations in the 
small avalanche region {siow < 20) to be observed. The 
peaks in the exponent correspond to even sizes. This 
indicates that small avalanches with odd sizes occur with 
a higher frequency than the frequency corresponding to 
a perfect power law. 

Note also that for the critical value Rc = 2.16 proposed 
in Ref. ,58,, the possible plateau at a height r' ^ 2.03 
does not exhibit a clear tendency to increase with L: 
it is very much distorted by odd-even fluctuations and, 
already for siow = 10^, has clearly decreased for all the 
studied system sizes. 

In our opinion, a final understanding of the behaviour 
of these distributions can only be achieved after a full 
finite-size scaling analysis of the distribution D{s;R,L) 
(involving the three variables), which is beyond the scope 
of this paper. Nevertheless, in view of the above observa- 
tions it seems that the distribution of avalanche sizes at 
the critical point shows two contributions: (i) the actual 
power-law behaviour corresponding to the non-spanning 
critical avalanches with t' ~ 1.7 that can only be ob- 
served at intermediate and large avalanche sizes, and (ii) 
the contribution from non-critical avalanches (lattice an- 
imals), in the small s region, which distorts the exponent 



towards the effective value r 



eff 



2.03. This scenario is 



not incompatible with previous numerical studies. 



VI. SUMMARY AND CONCLUSIONS 

The aim of this paper has been to illustrate the useful- 
ness of ML exponent maps in order to study distortions 
of power-law behaviour to the critical distributions of 
events. Such distortions are expected to occur for most 
experimental and numerical simulation data due to dif- 
ferent reasons. Fig. 1141 shows a schematic representation 
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of the main conclusions of this paper. On the left we show 
the ML exponent map that one can expect for experimen- 
tal data. When the sample set is small, a lack of statis- 
tics creates deformations of the theoretical plateau close 
to the diagonal axis of the map, which can be bounded 
by a proper estimation of the statistical error bars. Noise 
and undercounting of small-size events renders a defor- 
mation region starting in the bottom left-hand corner 
and which extends horizontally along the bottom border. 
Saturation of the amplifiers and counters also deforms 
the plateau in the upper right-hand corner and along the 
right edge. Between these three boundaries there should 
be a region with a well-defined plateau. The black dot 
indicates the best values of the high and low cut-off for 
the determination of the critical exponent. It should be 
mentioned that, if the size of the sample set is too small, 
such an ideal situation might not be achieved. 



Experimental Data Simulation Results 




FIG. 14. Schematic representation of ML exponent maps for 
experimental data (left) and simulation data (right). The 
lines indicate the regions where we can expect deformations 
of the theoretical plateau. The dot indicates the best values 
of the cut-offs for the estimation of the exponent. 



The right plot shows a similar scheme for numerical 
simulations of lattice systems. The deformation regions 
are similar to the previous ones. Along the bottom edge 
we find the footprint of the lattice character of the model 
that deforms the perfect power-law behaviour. Along the 
vertical right edge we find deformations associated with 
finite-size effects, which in some cases can be partially 
corrected. 



The examples analyzed in this work have included: 
three seismological catalogues corresponding to Japan, 
the San Andreas fault and El Hierro volcanic activity, 
three sets of Acoustic Emission data corresponding to the 
fracture of Vycor under compression, a cubic-tetragonal 
structural transition in FePd and a cubic-monoclinic 
structural transition in CuZnAl and, finally, numerical 
simulations of the 3D Random Field Ising model. These 
three case studies of the maps have allowed the good 
quality of the data to be checked or hypothesis for the 
observed distortions to be proposed. The studies have 
also been used to suggest improvements to measurements 
and/or numerical simulations. 
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